Published in : International Journal for Numerical Methods in Fluids (2006) 
Int. J. Numer. Meth. Fluids 2006; Vol 50: pp 421-436 



Fourth Order Compact Formulation of Navier-Stokes Equations 
and Driven Cavity Flow at High Reynolds Numbers 

E. Erturk* and C. Gokgol 

Energy Systems Engineering Department, Gebze Institute of Technology, 
^ Gebze, Kocaeli 414OO, Turkey 

o 

O SUMMARY 
D 

A new fourth order compact formulation for the steady 2-D incompressible Navier-Stokes equations is 
presented. The formulation is in the same form of the Navier-Stokes equations such that any numerical 
method that solve the Navier-Stokes equations can easily be applied to this fourth order compact 
formulation. In particular in this work the formulation is solved with an efficient numerical method that 
requires the solution of tridiagonal systems using a fine grid mesh of 601x601. Using this formulation, 
_i the steady 2-D incompressible flow in a driven cavity is solved up to Reynolds number of 20,000 with 

^ fourth order spatial accuracy. Detailed solutions are presented. 

KEY WORDS: High Order Compact Scheme; HOC; Steady 2-D Incompressible N-S Equations; PONS 
Equations; Driven Cavity Flow; High Reynolds Number Solutions 

o 

^ 1. INTRODUCTION 

High Order Compact (HOC) formulations are becoming more popular in Computational Fluid 
Dynamics (CFD) field of study. Compact formulations provide more accurate solutions in a 
^ compact stencil. 

O In finite differences, a standard three point discretization provides second order spatial 

j>. accuracy and this type of discretization is very widely used. When a high order spatial 

"'"] discretization is desired, ie. fourth order accuracy, then a five point discretization have to 

rN be used. However in a five point discretization there is a complexity in handling the points 

near the boundaries. 

High order compact schemes provide fourth order spatial accuracy in a 3 x 3 stencil and this 
type of compact formulations do not have the complexity near the boundaries that a standard 
wide (five point) fourth order formulation would have. 
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Dennis and Hudson [1], MacKinnon and Johnson [2], Gupta et al. [3], Spotz and Carey 
[4] and Li et al. [5] have demonstrated the efficiency of the high order compact schemes 
on the streamfunction and vorticity formulation of 2-D steady incompressible Navier-Stokes 
equations. 

In the literature, it is possible to find numerous different type of iterative numerical methods 
for the Navier-Stokes equations. These numerical methods, however, could not be easily used 
in HOC schemes because of the final form of the HOC formulations used in [1], [2], [3], [4] and 
[5]. This fact might be counted as a disadvantage of HOC formulations that the coding stage 
is rather complex due to the resulting stencil used in these studies. It would be very useful 
if any numerical method for the solution of Navier-Stokes equations described in books and 
papers could be easily applied to high order compact (HOC) formulations. 

In this study, we will present a new fourth order compact formulation. The difference of 
this formulation with References [1-5] is not in the way that the fourth order compact scheme 
is obtained. The main difference, however, is in the way that the final form of the equations 
are written. The main advantage of this formulation is that, any iterative numerical method 
used for Navier-Stokes equations, can be easily applied to this new HOC formulation, since 
the final form of the presented HOC formulation is in the same form with the Navier-Stokes 
equations. Moreover if someone already have a second order accurate {OAx'^) code for the 
solution of steady 2-D incompressible Navier-Stokes equations, using the presented formulation, 
they can easily convert their existing code to fourth order accuracy (OAa;^) by just adding 
some coefficients into their existing code. With this new compact formulation, we have solved 
the steady 2-D incompressible driven cavity flow at very high Reynolds numbers using a very 
fine grid mesh to demonstrate the efficiency of this new formulation. 



2. FOURTH ORDER COMPACT FORMULATION 

In non-dimensional form, steady 2-D incompressible Navier-Stokes equations in streamfunction 
{tp) and vorticity (w) formulation are given as 

1 d'^u _^ 1 d'^LO _di}dw _ dip du 
Re dx"^ Re dy'^ dy dx dx dy 

where x and y are the Cartesian coordinates and Re is the Reynolds number. For first order 
and second order derivatives the following discretizations are fourth order accurate 



d(j) _ 


^ Ax'^d^ 1 


dx 


6 dx^ 


d^c^ _ 


^ A.t2 d^4> 


dx'^ 


■^"^ 12 dx* ^ 



(3) 
(4) 



where (l)x and (^i^x are standard second order central discretizations such that 
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Pj+i - (Pi-i 
2Ax 

h+i - + 4>i-i 



If we apply the discretizations in Equations (3) and (4) to Equations (1) and (2), we obtain 
the following equations 

Ax'^d^^ Ay^aV 4, 
+ V'.. - ^ ^ - ^ ^ + 0{Ax\ A/) = -oj (7) 



1 1 1 Ax'^d^U 1 Ay'^d^UJ .^,,4,4, , 



Ay2 d^Tp Aa;2 , d^LO Ax'^ d^ip Ay'^ B^lo , 4 , , . 2 . 4^ .o^ 
-^'^^a^ - -^"^^9^ + ^'^^al + ^^^a^ + OiAx\Ax^Ay^Ay^) (8) 

In these equations we have third and fom-th derivatives {d^/dx^, d^/dy^, d'^/dx* and 
d'^/dy'^) of streamfunction and vorticity {tp and ui) variables. In order to find an expression for 
these derivatives we use Equations (1) and (2). For example, when we take the first and second 
x-derivative {d/dx and d"^ /dx^) of the streamfunction equation (1) we obtain the following 
equations 

d^ip _ duo d^tp 
dx^ dx dxdy"^ 

dx^ " ~dx^ ~ y > 

And also, by taking the first and second y-dcrivative {d/dy and d^/dy^) of the 
streamfunction equation (1) we obtain the following equations 

d^i) _ du d^tp ^^^^ 



dy^ dy dx'^dy 

d^ip d'^uj d^^p 



dy'^ dy"^ dx'^dy'^ 



(12) 



Using standard second order central discretizations given in Table I, these equations ((9), 
(10), (11) and (12)) can be written as the foUowings 
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dy- 



= -oJa: - i^xyy + 0(Aa;^ Ay^) (13) 

= -OJxx - ^xxyy + 0{Ax^,Ay^) (14) 

= -oJy - i'xxy + 0{Ax^, Ay'') (15) 

J = -^yy - i^xxyy + 0(Aa;2, Ay') (16) 



When we substitute Equations (14) and (16) into Equation (7) we obtain the following finite 
difference equation. 

Ax' Ay' ( Ax' Ay' \ 

^xx + ^yy = -0J- -J^^xx " -^^m " ( ^ + ^ j ^xxyy + 0{Ax*, Ax' Ay' , Ay*) (17) 

Wc note that, the sohition of Equation (17) is also a sohition to streamfunction equation 
(1) with fourth order spatial accuracy. Therefore if we numerically solve Equation (17), the 
solution we obtain will satisfy the streamfunction equation up to fourth order accuracy. 

In order to obtain a fourth order approximation for the vorticity equation (2), wc follow the 
same procedure. When we take the first and second derivatives of the vorticity equation (2) 
with respect to x- and ^-coordinates we obtain the followings 

d^u) ^ d'-tp doj ^ d^p d'oj ^ d'lp dui ^ d%p d'ui d^uj 
dx^ dxdy dx dy dx"^ dx' dy dx dxdy dxdy' 
d^oj _ d^tp duj d'lp d'uj d'i> d'oj j^^dip d^u> 
dx* dx'^dy dx dxdy dx"^ dxdy dx"^ dy dx^ 

^ d^ip du! ^ d'lp d'cu ^ d'pj d'uj ^ dtp d^uj d*uj 
dx^ dy dx"^ dxdy dx^ dxdy dx dx^dy dx'^dy'^ 

Qy?, Qy2 Qy, Qy Qxdy dxdy dy dx dy"^ dx'dy 



d*uj „ d^-tJj duj d'llj d'ljj d'lp d'u dip d^ 

Re-AT^+Re—^——+Re—-Pr-—-+Re 



UJ 



dy* dy^ dx dy^ dxdy dy'^ dxdy dy dxdy"^ 

^ d^tp dui ^ d'lp d'uj ^ d'lp d'oj ^ dtp d^uj d*uj 
dxdy' dy dxdy dy' dxdy dy' dx dy^ dx'dy' 

If we substitute Equations (18) and (20) for the third derivatives of vorticity {d'^oj/dx^ and 
d^uj/dy^) into Equations (8), (19) and (21) and also if wo substitute Equations (13) and (15) 
for the third derivatives of streamfunction {d^tp/dx^ and d^tp/dy^) into Equations (8), (19) 
and (21) and finally if we substitute Equations (19) and (21) for the fourth derivatives of 
vorticity {d*w/dx* and d*ui/dy*) into Equation (8), then we obtain the following equation 

Ax' Ay' Ax' Ay' 

^XX + ^yy — Re^-tPxyU>XX + Re—^tp^yiJyy + Re' —^tPytPyUJxX + RS.' tp ^tp ^LO yy = 
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' A.T^ Ay^ \ / A.x^ Ay^ 

RetpyUJ:^ - Rei^xUJy + Re ( —— + —— l/^xxyi^x - Re — — + — ^ ) Ipxyy^y 



Ax^ Aw^ Aa;^ Aw^ 

-Re^^^i'yi'xy^x + Re^ —^Ipx'tpyyUJx + Re^ —-tpytpxxi^y " Re^ ^^Ipxlpxy^y 



A.t2 Ay2\ /Aa;2 Ay^X Ax^ 

+ -J^ I ^y^xTO -Re{ — + — I ^xOJxxy - ^e — 



Ay2 ^^(l^x^ Ai/\ /Ax^ Ai/ 
- -^V'ra'^xj/ + -Re I — + — I i^xi'y'^xy - Re I— — 



+Re—;^^pyyLJxy + Re + "TFT I i^x'tpyi^xy - Re ( -p?; 7V ) WajWj^ 



12 12 



w,,^^ + 0(Ax*, Ax' Ay', Ay*) (22) 



Again wc note that, the solution of Equation (22) satisfy the vorticity equation (2) witli 
fourth order accuracy. 

As the final form of our HOC scheme, we prefer to write Equations (17) and (22) as the 
following 

i^XX +tpyy = -OJ + A (23) 

(1 + S) Uxx + (1 + ^vv = ^'^y + D)uJx- i^x + E)ivy + F (24) 

where 



— ^2 y2 I 12 12 ' ^^^yy 



Ax"^ Ax"^ 
B = -Re—^ipxy + Re"^ -^ipyipy 



2 



Ay2 Ay 

C = Re——i>xy + Re -^i'x'ip-. 



6 ^ 12 
Ax^ Aw2\ , ^ Aa;2 , , ^ Ay^ 
-j^ + ^ ) ^xxy - Re—tpyipxy + 1 

Aa;2 Ay^\ , „ Ax^ , , „ Ay^ 

+ ^ 1 ^xyy - Re — i'yipxx + — 1 



Ay2 /Aa;2 Ay^X /Aa;^ Ay^X 

+ + -Re —— + —— 1pxi^y(Jxy - -rz TTT '^^'^J/ 



6 " \ 12 12 y " " V 12 12 

"i^(ll" + l2")^^^^^ ^^^^ 



We note that, the finite difference Equations (23) and (24) are fourth order accurate 
(0(Aa;*, Aa;^Ay^, Ay*)) approximation of the streamfunction and vorticity equations (1) and 
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(2). In Equations (23) and (24), however, if A, B, C, D, E and F are chosen to be equal to 
zero then the finite difference Equations (23) and (24) simply become 

tpxx + i'yy = -W (26) 

-^Wxx + -^^yy = V'y'^x - V'xWy (27) 

Equations (26) and (27) arc the standard second order accurate (©(Aa;^, Ay^)) 
approximation of the streamfunction and vorticity equations (1) and (2). When we use 
Equations (23) and (24) for the numerical solution of 2-D steady incompressible Navier- 
Stokcs equations, wc can easily switch between second and fourth order accuracy just by 
using homogeneous values for the coefficients A, B, C, D, E and F or by using the expressions 
defined in Equation (25) in the code. 

In Equations (23), (24) and (25) instead of finite difference discretizations, if we substitute 
for partial derivatives we obtain the following differential equations 

^ _^ ^ _ _^ + ^ g 



Re dx^ Re dy"^ \dy J dx \dx J dy 

_Aa;2 92^ Ay^ d^cj / Ax^ Ay^\ 9^ 
rfda^ ^ 1^9j/2 " y-J^ dx^ 

B = -Re^^ + Re^^'^'t 
6 oxoy 12 dy dy 

6 dxdy 12 dx dx 

~ Vl^ 1^ ) dx'^dy ~ ^^dydxdy^ ^l^dxdy^ 
^ /Ax^ Ay'^\ d^ Ax'^ dj) d'^j) Ay^&ipdS^ 

\vr ^ IT ) dxdy"^ ^^Ihjlh^^ ^^dxdxdy 
Ax2 Ay^\d^l; d^uo /Ax^ Ay^\ dip d^uj Ax"^ d^ip d^u 



12 12 J dy dxdy^ V 12 12 / dx dx'^dy 6 dx"^ dxdy 

Ay^^d\j_ n Ay'^\d-<pdip d'^uj /Ax^ Ay'^\dujdw 

'^~Q~^ dxdy ^ ^' \ V2 ^ ~^ ) 'dx^ dxdy \W ~W ) dxdy 
1 /Ax2 Ay2\ ^4 



i?e V 12 12 y dx'^dy 

We note that the numerical solutions of Equations (28) and (29), strictly provided that 
second order discretizations in Table I are used and also strictly provided that a uniform 
grid mesh with Ax and Ay is used, are fourth order accurate to streamfunction and vorticity 
equations (1) and (2). We prefer to call Equations (28) and (29) Fourth Order Navier-Stokes 



HIGH ORDER COMPACT SCHEME 



7 



(FONS) equations. The only difference between FONS equations (28) and (29) and Navier- 
Stokes (NS) equations (1) and (2) are the coefficients A, B, C, D, E and F. In fact the NS 
equations are a subset of the FONS equations. We note that FONS equations (28) and (29) 
are in the same form with Navier-Stokes (NS) equations (1) and (2), therefore any iterative 
numerical method (such as SOR, ADI, factorization schemes, pseudo time iterations and etc.) 
used to solve streamfunction and vorticity equations (1) and (2) can also be easily applied to 
fourth order equations (28) and (29). Moreover, any existing code that solve the streamfunction 
and vorticity equations with second order accuracy can easily be modified to provide fourth 
order accuracy just by adding the coefficients A, B, C, D, E and F into the existing code to 
obtain the solution of FONS equations. Of course, when the coefficients A^ B, C, D, E and F 
are added into a second order accurate code to obtain fourth order accuracy, evaluating these 
coefficients would require extra CPU work. This might be considered as the cost of increasing 
accuracy from second order to fourth order. 



3. NUMERICAL METHOD 



Recently Erturk et al. [6] have presented a new, stable and efficient numerical method that solve 
the streamfunction and vorticity equations. The numerical method solve the governing steady 
equations through iterations in the pseudo time. In this study, we will apply the numerical 
method Erturk et al. [6] have proposed, to FONS equations (28) and (29) and solve the steady 
driven cavity flow with fourth order accuracy. For details about the numerical method, the 
reader is referred to Erturk et al. [6] . When we apply the numerical method to Equations (28) 
and (29), we obtain the following equations 



[l - A.^) (l - AtQ = r + At.- - AtA- + (a^^) (At^) r (31) 

1 - ^^^^ + iL ^ - + ^) " I ) = - 

E] ~~] (32) 
y -\ ' - ' jigQy2 ' " \^dx ' J dy / 

The solution methodology of these two equations are quite simple. First the streamfunction 
equation (31) is solved in two steps. For streamfunction equation, a new variable / is defined 
as the following 

1 - At^) = / (33) 





-) 


" d 




dx 


1 92 


+ At 


(dip 


' Redy'^ 


\ dx 



Using this variable in Equation (31) we obtain the following equation 
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1- Ai 



At 



(34) 



In this equation, the only unknown is the variable /. We first solve this equation for / by 
solving a tridiagonal system. After this, when we obtain the value of / at every grid point we 
solve Equation (33) for streamfunction (?/^"+^) by solving another tridiagonal system. 

After solving the streamfunction equation (31), we solve the vorticity equation (32). For 
this, similarly, we introduce a new variable g which is defined as the following 



1- At(l + C") 



1 92 



E 



d 



-At 

' Re dy^ \dx ' J dy 

Using this variable in Equation (32), we obtain the following equation 



(35) 



1 - At{l + B") 



1 



+ At 



+ At{l + B") 



dy 
1 



+ D 



9 = uj''- AtF" 



Atil + C" 



(36) 



d_ 

' Redx'^ ' ""\dy ' ^ ) dx ^ 

-At(^ dY — 
'Redx^ \dy J dx ^ 

) — —+At(^+EX — 
Re dy"^ \dx J dy ^ 

In this equation the only unknown is the variable g. By solving a tridiagonal system, we 
obtain the value of g ad, every grid point. Then we solve Equation (35) for vorticity (a;"+^) by 
solving another tridiagonal system. 

In a compact formulation, the stencil have 3x3 points. The solution at the first diagonal grid 
points near the corners of the cavity would require the vorticity values at the corner points. 
However, the corner points arc singular points for vorticity. Gupta et al. [7] have introduced an 
explicit asymptotic solution in the neighborhood of sharp corners. Similarly, Stortkuhl et al. [8] 
have presented an analytical asymptotic solutions near the corners of cavity and using finite 
clement bilinear shape functions they also have presented a singularity removed boundary 
condition for vorticity at the corner points as well as at the wall points. We follow Stortkuhl 
et al. [8] and use the following expression for calculating vorticity values at the wall 



3A/i2 



1 _4 1 
2^2 

111 



V 

T 



(37) 



where V is the speed of the wall which is equal to 1 for the moving top wall and equal to for 
the three stationary walls. For corner points, we use the following expression for calculating 
the vorticity values 



3A/i2 



• • 


• 






• • 


• 


• -2 


1 
2 

1 




1 

" 9 


• 1 


1 
2 
1 
4 



w = — 



V 
2h 



(38) 



where again V is equal to 1 for the upper two corners and it is equal to for the bottom two 
corners. The reader is referred to Stortkuhl et al. [8] for details. 
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4. RESULTS AND DISCUSSIONS 



The schematics of the driven cavity flow is given in Figure 1. In this figure the abbreviations 

BR, BL and TL refer to bottom right, bottom left and top left corners of the cavity, respectively. 
The number following these abbreviations refer to the vortices that appear in the flow, which 
are numbered according to size. 

For every Reynolds number considered, we have continued our iterations until, in the 
computational domain both the maximum residual of Equations (23) and (24), which are 
given as 



are less than 10~^°. Such a low value is chosen to ensure the accuracy of the solution. At these 
residual levels, the maximum absolute change in streamfunction value between two time steps, 
(max(|?/;"+-'^ — V-""!)); was in the order of lO"""^^ and for vorticity. (max(|a;"+-'^ — o^"!)), it was 
in the order of 10~^^. Obviously these convergence levels are far more less than satisfactory, 
however such low values demonstrate the efSciency of the numerical method used in this study 
which was presented by Erturk et al. [6]. 

Using an efficient numerical method, Erturk et al. [6] have clearly shown that numerical 
solutions of driven cavity fiow is computable for Re > 10,000 when a grid mesh larger than 
256x256 is used. With a grid mesh of 601x601 Erturk et al. [6] have solved the cavity flow 
up to i?e=21,000 using the numerical method also used in this study. In order to be able to 
obtain solutions at high Reynolds numbers, following Erturk et al. [6], in this study we have 
used a large grid mesh with 601 x 601 grids. With this many number of grid points we obtained 
steady solutions of the cavity flow up to i?e=20,000 with fourth order accuracy. 

Figures 2 to 6 show the streamfunction and vorticity contours of the driven cavity flow 
between i?e=l,000 and i?e=20,000. These figures show the vortices that are formed in the 
flow field as the Reynolds number increases. From these contour figures, we conclude that the 
fourth order compact formulation provides very smooth solutions. 

In Figure 7 we plot a very enlarged view of the top right corner (where the moving wall 
moves towards the stationary wall) of the streamfunction contour plot for the highest Reynolds 
number considered, i?e=20,000. In this figure the dotted lines show the grid lines. As it is seen 
in this enlarged figure, fourth order streamfunction contours are very smooth even at the first 
set of grid points near the corners. 

Table II tabulates the streamfunction and vorticity values at the center of the primary and 
secondary vortices and also the location of the center of these vortices for future references. 
This table is in good agreement with that of Erturk et al. [6] . 

Using Richardson extrapolation on the solutions obtained with different grid meshes, 
Erturk et al. [6] have presented theoretically fourth and sixth order accurate (OAa;^ and 




(39) 





(40) 
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OAx^) strcamfunction and vorticity values at the center of the primary vortex. Table III 
and IV compares the forth order compact scheme solutions of the streamfunction and the 
vorticity values at the center of the primary vortex with the fourth order (OAa;^) Richardson 
extrapolated solutions tabulated in Erturk et al. [6]. The present solutions and the solutions 
of Erturk et al. [6] agree with each other. 

5. CONCLUSIONS 

In this study a new fourth order compact formulation is presented. The uniqueness of this 
formulation is that the final form of the HOC formulation is in the same form of the Navier- 
Stokes equations such that any numerical method that solve the Navier-Stokes equations can 
be easily applied to the PONS equations in order to obtain fourth order accurate solutions 
[O/^x'^). Moreover with this formulation, any existing code that solve the Navier-Stokes 
equations with second order accuracy (OAa;^) can be altered to provide fourth order accurate 
(OA.x'*) solutions just by adding some coefficients into the code at the expense of extra CPU 
work of evaluating these coefficients. 

In this study, the presented fourth order compact formulation is solved with a very efficient 
numerical method introduced by Erturk et al. [6]. Using a fine grid mesh of 601x601, as it 
was suggested by Erturk et al. [6] in order to be able to compute for high Reynolds numbers, 
the driven cavity flow is solved up to Reynolds numbers of i?e=20,000. The solutions obtained 
agree well with previous studies. The presented fourth order accurate compact formulation is 
proved to be very efficient. 
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Figure 1. Schematic view of driven cavity flow 



Figure 2. Streamfunction and vorticity contours for Re=1000 



Figure 3. Streamfunction and vorticity contours for Re=5000 



Streamfunction and vorticity contours for Re=10000 




Figure 5. Streamfunction and vorticity contours for Re=15000 




Figure 6. Streamfunction and vorticity contours for Re=20000 




Figure 7. Streamfunction contours for Re=20000, enlarged view of top right comer. 
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Table 1. Standard Second Order Central Discretizations, 0{Ax'^ , Ay"^). 



Re 

1000 2500 5000 7500 10000 12500 15000 17500 20000 



Primary tl? 
Vortex oj 

> y) 


-0.118938 
-2.067760 
(0.5300 , 0.5650) 


-0.121472 
-1.976132 
(0.5200 , 0.5433) 


-0.122216 
-1.940547 
(0.5150 , 0.5350) 


-0.122344 
-1.926478 
(0.5133 , 0.5317) 


-0.122306 
-1.918187 
(0.5117 , 0.5300) 


-0.122201 
-1.912307 
(0.5117 , 0.5283) 


-0.122060 
-1.907651 
(0.5100 , 0.5283) 


-0.121889 
-1.903659 
(0.5100 , 0.5283) 


-0.121694 
-1.900032 
(0.5100 , 0.5267) 


BRl u 

> y) 


0.17297B-02 
1.118222 
(0.8633 , 0.1117) 


0.26623E-02 
1.955594 
(0.8333 , 0.0900) 


0.30735E-02 
2.739071 
(0.8050 , 0.0733) 


0.32265E-02 
3.243925 
(0.7900 , 0.0650) 


0.31896B-02 
3.756425 
(0.7750 , 0.0600) 


0.30972E-02 
4.357323 
(0.7600 , 0.0550) 


0.30022B-02 
4.965304 
(0.7450 , 0.0500) 


0.29021E-02 
5.568522 
(0.7333 , 0.0467) 


0.28012B-02 
6.125275 
(0.7200 , 0.0433) 


BLl u> 

. y) 


0.23345B-03 

0.354271 
(0.0833 , 0.0783) 


0.93093E-03 

0.979831 
(0.0833 , 0.1117) 


0.13758B-02 

1.514292 
(0.0733 , 0.1367) 


0.15337E-02 

1.858636 
(0.0650 , 0.1517) 


0.16135B-02 

2.182043 
(0.0583 , 0.1633) 


0.16568E-02 

2.358513 
(0.0550 , 0.1683) 


0.16663B-02 

2.490168 
(0.0533 , 0.1717) 


0.16450E-02 

2.704118 
(0.0500 , 0.1767) 


0.16083B-02 

2.885054 
(0.0483 , 0.1817) 


BR2 u 

(5<; . y) 


-0.49242E-07 
-0.76887B-02 
(0.9917 , 0.0067) 


-0.12035E-06 
-0.1193SE-01 
(0.9900 , 0.0083) 


-0.14271E-05 
-0.33647E-01 
(0.9783 , 0.0183) 


-0.32742E-04 
-0.155445 
(0.9517 , 0.0417) 


-0.14014E-03 
-0.309618 
(0.9350 , 0.0683) 


-0.25498E-03 
-0.398536 
(0.9283 , 0.0817) 


-0.34006E-03 
-0.456877 
(0.9267 , 0.0883) 


-0.40492E-03 
-0.512862 
(0.9283 , 0.0967) 


-0.46117E-03 
-0.554507 
(0.9300 , 0.1050) 


BL2 u 

('■^ . y) 


-0.65156E-08 
-0.29579E-02 
(0.0050 , 0.0050) 


-0.26763E-07 
-0.93547E-02 
(0.0067 , 0.0067) 


-0.65524E-07 
-0.12183E-01 
(0.0083 , 0.0083) 


-0.20189E-06 
-0.16587B-01 
(0.0117 , 0.0117) 


-0.11017E-05 
-0.30153E-01 
(0.0167 , 0.0200) 


-0.64493E-05 
-0.74310B-01 
(0.0267 , 0.0317) 


-0.22287E-04 
-0.140685 
(0.0383 , 0.0417) 


-0.48296E-04 
-0.200227 
(0.0500 , 0.0483) 


-0.78077E-04 
-0.243259 
(0.0583 , 0.0533) 


TLl w 

('■^ . y) 




0.34284E-03 
1.344040 
(0.0433 , 0.8900) 


0.14459B-02 
2.115544 
(0.0633 , 0.9100) 


0.21311E-02 
2.236986 
(0.0667 , 0.9117) 


0.26237B-02 
2.322310 
(0.0700 , 0.9100) 


0.29949E-02 
2.375671 
(0.0733 , 0.9100) 


0.32865B-02 
2.424136 
(0.0767 , 0.9100) 


0.35245E-02 
2.463058 
(0.0783 , 0.9117) 


0.37247B-02 
2.497615 
(0.0800 , 0.9117) 


BR3 u 

(x > y) 


_ 




0.42053B-10 

0.66735E-03 
(0.9983 , 0.0017) 


0.75831B-09 

0.70882E-03 
(0.9967 , 0.0017) 


0.37451B-08 

0.24926E-02 
(0.9967 , 0.0050) 


0.78212E-08 

0.37843E-02 
(0.9950 , 0.0050) 


0.11742B-07 

0.33420E-02 
(0.9950 , 0.0050) 


0.16951E-07 

0.42443E-02 
(0.9950 , 0.0067) 


0.27181B-07 

0.49718E-02 
(0.9933 , 0.0067) 


BL3 u 

(5<; ; y) 


- 


_ 








0.24668E-09 
0.50277E-03 

(U.UU17 , U.UU17J 


0.57300E-09 
0.10067E-02 

(^U.UUl? , U.UUoo; 


0.13069E-08 
0.21872E-02 

(^U.UUOO , U.UU33J 


0.23598E-08 
0.18079E-02 

tU.UU33 , U.UUooJ 


TL2 w 

(x , y) 












-0.16253B-05 
-0.204044 
(0.0067 , 0.8300) 


-0.15848E-04 
-0.524342 
(0.0150 , 0.8250) 


-0.41183B-04 
-0.737721 
(0.0200 , 0.8217) 


-0.70762E-04 
-0.982968 
(0.0250 , 0.8200) 



Table 2. Properties of primary and secondary vortices; streamfunction and vorticity values, (x,y) locations 



Re 


present 


Erturk et. al. [2] 


1000 


-0.118938 


-0.118939 


2500 


-0.121472 


-0.121469 


5000 


-0.122216 


-0.122213 


7500 


-0.122344 


-0.122341 


10000 


-0.122306 


-0.122313 


12500 


-0.122201 


-0.122229 


15000 


-0.122060 


-0.122124 


17500 


-0.121889 


-0.122016 


20000 


-0.121694 


-0.121901 



Table 3. Minimum Streamfunction values at the Primary Vortex 
for various Reynolds numbers 



Re 


present 


Erturk et. al. [2] 


1000 


-2.067760 


-2.067579 


2500 


-1.976132 


-1.976096 


5000 


-1.940547 


-1.940451 


7500 


-1.926478 


-1.926282 


10000 


-1.918187 


-1.917919 


12500 


-1.912307 


-1.912072 


15000 


-1.907651 


-1.907602 


17500 


-1.903659 


-1.903975 


20000 


-1.900032 


-1.900891 



Table 4. Vorticity values at the center of the Primary Vortex 
for various Reynolds numbers 



